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Low-mass protostars form from condensations inside molecular clouds when gravity 
overwhelms thermal and magnetic supporting forces. The first phases of the formation of a 
solar-type star are characterized by dramatic changes not only in the physical structure but 
also in the chemical composition. Since PPIV (e.g., Langer et al.), exciting new developments 
have occurred in our understanding of the processes driving this chemical evolution. These 
developments include two new discoveries : 1) extremely enhanced molecular deuteration, 
which is caused by the freeze-out of heavy-element-bearing molecules onto grain mantles 
during the Prestellar Core and Class source phases; and 2) hot corinos, which are warm 
and dense regions at the center of Class source envelopes and which are characterized by a 
multitude of complex organic molecules. In this chapter we will review these two new topics, 
and will show how they contribute to our understanding of the first phases of solar-type stars. 
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1. INTRODUCTION 

Molecular clouds are the placentas inside which matter 
evolves from embryos to stars and planetary systems. The 
material in the molecular cloud feeds the newly forming 
system. During this formation, the material undergoes sev- 
eral dramatic changes. This chapter focuses on the changes 
of the chemical composition during the first phases of star 
formation. 

There are several reasons why the study of chemistry in 
the first phases of star formation is fascinating and impor- 
tant. Among them, two stand out. First, chemistry is a 
very powerful diagnostic, both of the current and the past 
physical conditions of the forming protostar. Much of this 
chapter will be spent in illustrating that point in detail. In 
addition, chemistry in the first phases of star formation may 
affect the chemical composition of the objects that will even- 
tually form the planetary system: planets, comets, and as- 
teroids. In this sense, the study of chemistry during the first 
phases of star formation is far-reaching. Indeed, one of the 
major and more fascinating questions linked to the process 
of the formation of a star and its planetary system, espe- 
cially if similar to our Solar System, is: what is the chemical 



budget acquired during the protostellar phase and inher- 
ited by the forming planets? Answering this question might 
even shed some light on our own origins. The above ulti- 
mate question implies answering several linked questions: 
what molecules are formed during the protostellar phase? 
Do molecules exist that are formed prevalently during the 
protostellar phase, and that are therefore a hallmark of this 
period? What is their fate? Do they condense onto the grain 
mantles during the proto-planetary phase? Are they incor- 
porated into planetesimals, which eventually form comets, 
meteorites and planets? Are these pristine molecules re- 
leased into nascent planetary atmospheres during early, in- 
tense cometary bombardment? And what is the ultimate 
molecular complexity reached during star formation? In 
this contribution we will focus on the first questions, those 
related to the first phase of star formation. 

The story starts with the first step towards collapse: the 
Prestellar Core phase (see the chapters by Di Francesco et 
al. and Ward-Thompson et al). These objects are cold 
(< 10 K) and dense (> 10 5 cm~ 3 ) condensations inside 
molecular clouds. What makes them particularly impor- 
tant is that they are believed to be on the verge of col- 
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Table 1 : Summary of the properties of a sample of Prestellar 
Cores. n c (H 2 ) is the central density; / D (CO) is the CO 
depletion factor; the Molecular D/H ratio column contains 
the observed N 2 D + /N 2 H+ and/or (D2CO/H2CO) 1 / 2 ratio. 



Name 


n c (H 2 ) 


/d(CO) 


Mol. 


Ref. 




(10 5 cm" 3 ) 




D/H 


a 


L1521F 


3 


2.5 


0.01 


1,2,3 


B68 


1 


3.4 


0.03 


3,4 


L1689B 


1 


4.5 


0.09 


3,6,7,8 


L183 


10 


12 


0.22 


3,9 


LI 544 


14 


14 


0.23 


3,6,10 


L694-2 


9 


11 


0.26 


3 


L429 


6 


16 


0.28 


3,6,7 


OphD 


3 


14 


0.44 


3,6,7 


LI 709 A 


1 


5.5 


0.17 


6,7 


L1498 


1 


7.5 


0.04 


3,5 


L1517B 


2 


9.5 


0.06 


3,5 



"REFERENCES: 1- Tafalla and Santiago (2004); 2- Hirota 
et al. (2001); 3- Crapsi et al. (2005a); 4- Bergin et al. 
(2001); 5- Tafalla et al. (2004); 6- Bacmann et al. (2002); 
7- Bacmann et al. (2003); 8- Redman et al. (2002); 9- 
Pagani et al. (2005); 10- Tafalla et al. (2002). 

lapse (Tafalla et al., 1998; Caselli et al, 2002b; Young 
et al, 2004; Crapsi et al, 2005b), and, in this sense, they 
are considered to be representative of the initial conditions 
of star formation. During this phase, matter slowly accu- 
mulates towards the center under the gravitational force, 
which counteracts the thermal and/or magnetic pressure. 
As the density increases, gaseous molecules start to freeze- 
out onto the cold dust grains, forming H20-dominated ice 
mantles "dirtied" with several other molecules. The pro- 
cess is so efficient that in the innermost and densest regions 
of the condensation, heavy-element-bearing molecules are 
thought to be virtually all frozen out. The low temperatures 
and the disappearance of most molecules, and particularly 
of CO, from the gas phase trigger a peculiar chemistry: an 
extreme molecular deuteration. 

Once gravitational contraction finally takes over, a proto- 
star is born, consisting of a central object, which eventually 
will become a star, surrounded by an envelope from which 
the future star accretes matter. In the beginning, the enve- 
lope is so thick that it obscures the central object, and the 
Spectral Energy Distribution (SED) is totally dominated by 
the cold outer regions of the envelope, with a temperature 
lower than ~ 30 K. This phase is represented by so-called 
Class sources (Andre et al, 2000). Most of the envelope is 
cold and depleted of heavy-bearing molecules, frozen onto 
the grain mantles, exactly as in Prestellar Cores. However, 
the presence of a central source, powered by gravitational 
energy, causes the heating of the innermost regions of the 
envelope. In these regions the dust temperature reaches 100 
K, causing the evaporation of the grain mantles formed dur- 
ing the Prestellar-Core phase. The molecules trapped in the 
ices are injected back into the gas phase, giving rise to a 



rich and peculiar chemistry. Complex organic molecules 
are found in these regions, called Hot Corinos. 

In this chapter, we will show that extreme deuteration 
and Hot Corinos are specific signatures of the first phases 
of low-mass star formation. The chapter is organized as 
follows. We first briefly describe the physical and chemical 
structures of Prestellar Cores and Class source (Section 
2). We then discuss in detail the deuteration phenomenon 
(Section 3), Hot Corinos (Section 4), and the chemical mod- 
els that have been developed to explain these observations 
(Section 5). The chapter ends with our conclusions. 

We conclude this section by listing a few reviews ap- 
peared in the literature related to the subjects treated in this 
chapter. Previous reviews of the chemistry in general are 
reported in Herbst (2005), and during the star formation pe- 
riod in van Dishoeck and Blake (1998), Longer et al. (2000) 
and Caselli (2006). Reviews of the observations of mul- 
tiply deuterated molecules as well as the related chemical 
models can be found in Roueff and Gerin (2003), Cecca- 
relli (2004a) and Millar (2005), and on the Hot Corinos in 
Ceccarelli (2004b). 

2. THE PHYSICAL AND CHEMICAL STRUC- 
TURE IN THE FIRST PHASES OF THE COL- 
LAPSE 

Since PPIV, our understanding of the physical (temper- 
ature and density profiles) and chemical (molecular and 
atomic abundance profiles) structure of matter in the first 
phases of low-mass star formation has improved consider- 
ably. This is due both to the increased sensitivity of ground- 
based instruments and to the development of sophisticated 
models to interpret the data. In this section, we review 
the observations and their interpretations, which have led 
to the reconstruction of the gas density, temperature and 
molecular composition in both Prestellar Cores and Class 
sources. 

2.1. Prestellar Cores 

Several authors have studied the physical and chemical 
structure of starless cores; detailed reviews are presented in 
the chapters by Di Francesco et al. and Ward-Thompson et 
ah. Here we briefly recall the basic properties characteriz- 
ing Prestellar Cores, and discuss some aspects complemen- 
tary to those presented by the mentioned authors. First, the 
density profiles of the studied Prestellar Cores are reason- 
ably well represented by Bonnor-Ebert spheres. In practice, 
the density profile can be approximated by a power law r~ 2 
in the outer regions and a plateau in the center regions of 
the Prestellar Core. In the following we will refer to the 
density in the plateau as the "central" density. Typically, 
the radius of the plateau is 3000-6000 AU (Crapsi et al, 
2005a). Figure 1 shows the example of LI 544, one of the 
best studied Prestellar Cores (Tafalla et al, 1998; Caselli 
et al, 2002b, c), but other Prestellar Cores present basically 
the same structure (Alves et al, 2001; Crapsi et al, 2005a). 
In the same figure, the predicted temperature profile is also 
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Fig. 1. — The density and temperature profiles derived for 
the Prestellar Core LI 544 (adapted from Crapsi et al. in 
prep.). 

shown. Note the temperature drop in the inner 3000 AU 
region, confirmed by observations of NH 3 {Craspi et al. 
in prep). A similar drop in the dust temperature has been 
previously observed in other Prestellar Cores {Evans et al, 
2001; Pagani et al, 2003, 2004, 2005), as expected from 
radiative transfer modeling of centrally concentrated and 
externally heated dense cloud cores (Zucconi et al., 2001). 
These large central densities and low temperatures are ac- 
companied by the depletion of CO molecules in the central 
regions, as they condense out onto the grain mantles (Evans 
et al, 2001; Caselli et al, 1999; Galli et al, 2002). Indeed, 
if one had to give a short list of properties for the Prestellar 
Cores they would be large central densities: > 10 5 cm -3 ; 
very low central temperatures: < 10 K; depletion of molec- 
ular species, including CO; and enhanced molecular deuter- 
ation. Table 1 summarizes these properties in a sample of 
Prestellar Cores. In Figure 2 we show the observed CO 
depletion factors versus the central densities, where the de- 
pletion factor is computed with respect to the standard CO 
abundance in molecular clouds, 9.5 x 10~ 5 with respect to 
H 2 (Frerking et al, 1982). As can easily be seen, there is a 
clear correlation between the CO depletion factors and the 
densities. This observationally confirms the basic idea that 
CO molecules disappear from the gas phase because they 
freeze-out onto the grains, and that the condensation rate 
is proportional to the density (see below). Figure 2 also 
reports theoretical curves, obtained assuming that the CO 
molecules condense onto the grains at a rate fcf ree ze and are 
released in the gas phase because of the evaporation caused 
by cosmic rays at a rate fc cr = 9.8 x 10~ 15 s -1 (Hasegawa 
andHerbst, 1993). The freezing rate can be written as fol- 
lows: 

fcfrcozo = S (na^vco) n g , (1) 

where we adopted a sticking coefficient S = 1 (Burke 
and Hollenbach, 1983), and a mean grain radius of a gI . 



o 
o 



c 
o 



10 
_cu 

CL 
CD 

^ 5 

o 
o 







10 


i 1 o 5 / no 4 ; 






/ * / 






/ * / * 
1 J 






1 M 






* / / : 




/* 






/ * 






/ * 













1 10 
Density ( 1 5 cm~ 3 ) 

Fig. 2. — CO depletion factor versus the H 2 density. The 
asterisks represent the observations in the Prestellar Cores 
quoted in Table 1. The solid lines show the theoretical 
curves for an age of 1 x 10 4 and 1 x 10 5 yr, respectively, as- 
suming an average grain radius of 0.1 fim. The dashed lines 
are computed taking an average grain radius of 0.5 fim. 

The grain number density n g is given by the (mass) dust-to- 
gas ratio (0.01) multiplied by the gas density, and divided 
by the grain mass (computed assuming a grain density of 
2.5 gm cm -3 ). The figure shows curves assuming an av- 
erage grain radius of 0. 1 (solid lines) and 0.5 p,m (dashed 
lines) respectively. In the first case, which is the typical 
value assumed in chemical models for the ISM, the denser 
Prestellar Cores lie around the curve with an age of 1 x 10 4 
yr, which is definitively too short an age for these objects, 
based on several arguments. However, if the average grain 
radius is larger and equal to 0.5 /im, for example, the ob- 
served points lie around the curve at 1 x 10 5 yr, which is a 
more realistic estimate of the age of these objects. Larger 
grains would shift the curve even more to larger ages. Note 
that we used the same value for the cosmic ray evaporation 
rate, although increasing the grain sizes decreases it because 
of the larger volume to heat (Bringa and Johnson, 2004; 
Shen et al, 2004), enhancing the effect. Since the surface 
7rcig r is dominated by the small grains, Figure 2 suggests 
that the small grains are efficiently removed by coagula- 
tion with large grains in the innermost regions of Prestellar 
Cores (Ossenkopf, 1993; Bianchi et al, 2003; Flower et al, 
2005). Other explanations are also possible. For example, 
the cosmic ray evaporation rate can be enhanced if the ice 
mantles contain a substantial ( 1 %) fraction of radicals due 
to UV photolysis (D'Hendecourt et al, 1982, 1986). How- 
ever, in order to explain Figure 2 the mantle evaporation 
rate should be increased by more than a factor 30 with re- 
spect to the Hasegawa and Herbst (1993) rate. Therefore, 
our preliminary conclusion is that, most likely, the average 
sizes of the grains in Prestellar Cores are larger than in the 
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ISM, and that these objects are older than about 1 x 10 5 yr. 
It is also interesting noticing that less dense Prestellar Cores 
point to greater ages and/or smaller average grain sizes. It 
is not clear at this stage whether this means that those cores 
will not evolve into denser ones, or whether Prestellar Cores 
spend more time in this first "low-density" phase than in the 
"high-density" phase. 

2.2. Class sources 

The physical and chemical structure of the envelopes of 
Class sources is affected by the presence of a central ob- 
ject, i.e. a heating source, which implies gradients in the 
density and temperature. This, in turn, implies a gradi- 
ent in the chemical composition of the gas across the en- 
velope. Two basic methods have been used in the litera- 
ture to derive the chemical and physical structure of Class 
envelopes: i) modeling of multi-frequency observations of 
several molecules, and ii) modeling of the continuum SED 
and map, coupled with molecular multi-frequency obser- 
vations. The first method makes use of the intrinsic prop- 
erty that each line from any molecule probes the specific 
region where the line is excited, namely a region with a 
specific density and temperature. Using several lines from 
the same molecule can, therefore, be used to reconstruct 
the gas density and temperature profiles of the envelope. 
The observations of several molecules are then used to de- 
rive the abundance profile of each molecule. Theoretical 
examples of this method, applied to the problem of collaps- 
ing envelopes, are discussed in Ceccarelli et al. (2000a), 
Maret et al. (2002a), Ceccarelli et al. (2003), and Parise 
et al. (2005b). In these models the density is assumed to 
follow the evolution of a Singular Isothermal Sphere con- 
tracting under the hypothesis of isothermal collapse - the 
so-called "inside-out" framework, developed by Shu and 
collaborators {Shu, 1977). In this framework, the density 
follows a two-slope power law: in the regions not reached 
by the collapse the density has an r -2 dependence, whereas 
in the collapsing regions the density follows an r~ 3 / 2 law. 
In addition, in using molecular line observations, one has to 
also take into account the abundance profile of the molecule 
used. The simplest profile is that of a single step function, 
which assumes that the relevant molecule has an almost 
constant abundance in the cold outer envelope and a jump 
in the inner region when the grain mantles sublimate. Other 
profiles have been used though, in which the abundance also 
has a drop in the outer envelope, to take into account the 
molecular depletion as in Prestellar Cores {Ceccarelli et al, 
2001; Schdier et al, 2004; Jdrgensen et al., 2004). 

The second method uses the dust continuum SED si- 
multaneously with maps to trace the H2 distribution, and 
the observations are best-fitted with assumed single power 
law density distributions {Shirley et al, 2002; J0rgensen 
et al., 2002; Williams et al., 2005; Young and Evans, 2005). 
Then, in order to derive the molecular abundance profiles, 
single-jump and/or drop models are used to interpret multi- 
frequency observations, as in Schdier et al. (2002), Schdier 
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Fig. 3. — The temperature (upper panel) and density (lower 
panel) profile of the envelope surrounding the Class 
protostar IRAS 16293-2422 (figure from Ceccarelli et al. 
(2000a)). 



et al. (2004), Jprgensen et al. (2004), and Evans et al. 
(2005). Both methods give approximatively the same struc- 
ture (temperature and density distribution) when applied to 
the same object, as in the case of IRAS 16293-2422 {Cecca- 
relli et al, 2000a; Schdier et al, 2002). Figure 3 shows the 
derived structure of the envelope of this source, obtained by 
applying the first method to several molecules {Ceccarelli 
et aZ.,2000a,b). 

Both classes of models predict the existence of regions 
with dust temperatures larger than 100 K, the sublimation 
temperature of the grain mantles {Ceccarelli et al, 1996). 
In these regions, therefore, the components of the grain 
mantles are injected into the gas phase, and the abundances 
of the corresponding molecules should increase with re- 
spect to the abundance in the cold {Tdust < 100 K) enve- 
lope. Water is the major component of the ice mantles, but 
unfortunately water lines are not observable from ground- 
based telescopes. The spectrometer LWS on board ISO 
provided observations of water lines towards a few Class 
sources {Ceccarelli et al, 1999), but with relatively poor 
spatial and spectral resolution, so that the interpretation of 
the data is not unique. Two interpretations have been ad- 
vanced. The first one assumes that the water lines orig- 
inate in the shocks at the interface between the outflows 
and the envelopes of these sources {Giannini et al., 2001; 
Nisini et al., 2002). The second one assumes that the wa- 
ter lines originate mostly in the envelope. Based only on 
the ISO data, it is impossible to discriminate observation- 
ally between the two hypotheses. However, if the water line 
spectra are interpreted as being emitted in the envelope, al- 
lowing a jump in the water abundance in the Tdust > 100 K 
region, the agreement among the model predictions and the 
observed line fluxes is rather good {Ceccarelli et al., 2000a; 
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Maret et al, 2002b). The observed data are reproduced if 
the water abundance is ~ 3 x 10~ 7 in the outer envelope, 
and jumps by about a factor 10 in the inner region. In sup- 
port of this interpretation, the physical structure derived by 
this analysis is substantially confirmed by observations of 
other molecules obtained with ground-based telescopes (see 
below). 

Formaldehyde is also an important grain mantle compo- 
nent, and this molecule has the advantage that it has several 
transitions in the millimeter to sub-millimeter wavelength 
range observable with ground-based telescopes (Mangum 
and Wootten, 1993; Ceccarelli et al, 2003). These observa- 
tions have a much better spatial and spectral resolution than 
ISO-LWS and provide stronger constraints for the models. 
A survey of the H 2 CO line emission towards a bit fewer 
than a dozen Class sources has shown that all the tar- 
geted sources, except VLA16293, have a region where the 
formaldehyde abundance jumps by more than one order of 
magnitude (Ceccarelli et ai, 2000b; Maret et ai, 2004). 
Several of these sources also show jumps in the methanol 
abundance (Maret et ai, 2005). The predicted sizes of the 
jump regions range from 10 to 150 AU, and are, therefore, 
comparable to the sizes of the Solar System. The densi- 
ties are predicted to be larger than 10 8 cm~ 3 . However, 
one should be aware that the above observations have been 
obtained with single-dish telescopes, which, at best, have 
spatial resolutions corresponding to about 1000 AU in ra- 
dius. Therefore, both the estimates of the abundances and 
the warm region sizes suffer from a relatively large uncer- 
tainty. Besides, the details of the model adopted in the 
analysis of the data give rise to additional systematic un- 
certainties so that the very existence of the H2CO abun- 
dance jumps is contested by some authors (Schdier et ai, 
2004; J0rgensen et ai, 2005c, a). In addition, the situation 
is further confused by the presence at small scales of cavi- 
ties and outflows, which some authors also consider a ma- 
jor component in the line emission attributed to the abun- 
dance jumps of the above analysis (Chandler et ai, 2005; 
Jfirgensen et ai, 2005b). 

Finally, several simple molecules ( HCO+, N 2 H + , CS, 
SO, S0 2 , HCN, HNC, HC 3 N and CN) have been observed 
in a sample of Class sources (Jfirgensen et ai, 2004). 
Given the observed transitions, this study was able to trace 
the abundance of the observed species in the outer enve- 
lope and probe some formation routes, such as for N 2 H + 
and HCO + . One remarkable result is the discovery of large 
regions in the outer envelopes where the molecular abun- 
dances are depleted, because of condensation onto grain 
mantles (Maret et ai, 2004; J0rgensen et al., 2005d). How- 
ever, even more remarkable is the discovery of extreme 
deuterium fractionation in the Class sources, which will 
be discussed in Section 3, and which is indeed linked to the 
molecular depletion. 

In summary, from a chemical point of view, the similar- 
ity to Prestellar Cores is the most prominent characteristic 
of the outer envelopes of Class sources. With respect to 



molecular depletion and extreme deuteration, they are in- 
deed virtually indistinguishable, suggesting that Prestellar 
Cores are the likely precursors of Class sources (Maret 
et ai, 2004; J0rgensen et al., 2005c). 

3. EXTREME DEUTERATION 

The previous section has anticipated a major character- 
istic shared by Prestellar Cores and the envelopes of Class 
sources: an extreme molecular deuteration. Although 
the deuterium abundance is 1.5 x 10~ 5 relative to hydro- 
gen (Linsky, 2003), singly deuterated molecules have been 
observed in both types of objects with abundances relative 
to their hydrogenated counterparts between 10% and 50% 
(Guelin et al, 1977; Loren et al, 1990; van Dishoeck et al, 
1995; Tine et al, 2000; Shah and Wootten, 2001; Roberts 
et al, 2002; Crapsi et al, 2004, 2005a). Even more ex- 
treme, doubly and triply deuterated molecules have been 
observed (Table 2) with D/H ratios reaching 30% for D 2 CO 
(Loinard et al, 2002) and 3% for CD3OH (Parise et al, 

2004) . This enhances the molecular D/H ratio by up to 13 
orders of magnitude with respect to the elemental D/H ra- 
tio. In this section, we review these spectacular observa- 
tions, which have led to the development of a new class of 
models for molecular deuteration, described in Section 5. It 
is worth emphasizing that extreme fractionations have only 
been observed in low mass Prestellar Cores (Section 3.1) 
and protostars (Section 3.2) so far. Massive protostars do 
not show the same deuterium enrichment, probably because 
the chemistry occurs in warmer environments. 

3.1. Prestellar Cores 

In the last few years, it has become increasingly clear 
that molecular deuteration reaches extreme values in 
Prestellar Cores, where multiply deuterated molecules have 
been detected with ratios larger than 1%. The list of multi- 
ply deuterated molecules detected in Prestellar Cores (and 
Class sources) is reported in Table 2, as well as the typi- 
cal measured D/H ratio for each species. To the best of our 
knowledge, two surveys of multiply deuterated molecules 
- D 2 CO (Bacmann et al, 2003) and ND 3 (Roueff et al, 

2005) - have been obtained in (small) samples of Prestellar 
Cores. Bacmann et al. (2003) found that the D2CO/H2CO 
ratios correlate with the CO depletion factors, proving ob- 
servationally the key role of CO depletion in the deuterium 
enhanced fractionation. This trend is confirmed now by the 
larger sample of Table 2, as shown in Figure 4: the larger 
the CO depletion the larger the molecular D/H ratio (see 
also Crapsi et al. (2005a)). This is a key point whose impor- 
tance will be further discussed in Section 5. Another aspect 
shown in the figure is the similarity of the D/H ratio of the 
three observed molecules: N 2 H + , H 2 CO and NH3. Note 
that the three molecules refer to the D/H ratios of singly, 
doubly and triply deuterated isotopologues respectively, 
and the ratios shown in the figure are scaled to account for 
the number of D atoms. While N 2 H+ is certainly a gas 
phase product (formed by the reaction of N 2 with H3 ), 
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Table 2: Summary of the multiply deuterated molecules ob- 
served in Prestellar Cores and Class sources. Typical deu- 
terium fractionation ratios are given with respect to the fully 
hydrogenated counterparts, except for D2S and HD \ where 
the ratio is with respect to the singly deuterated isotopo- 
logue. Note that variations among different sources amount 
to about a factor 3 or so. 



Species 


Deuterium 


Ref. 




Fractionation 


a 


NHD 2 


0.03 


1,2,3 


ND 3 


0.006 


3,4,5 


D 2 CO 


0.15 


6,7,8 


D 2 S 


0.12 


9 


HD+ 


1 


10 


CHD 2 OH 


0.06 


11,8 


CD3OH 


0.03 


12 


D 2 CS 


0.10 


13 



"REFERENCES: 1- Roueff et al (2000); 2- Loinard et al. 
(2001); 3- Roueffetal. (2005) 4- Lis et al. (2002); 5- van der 
Tak et al. (2002); 6- Ceccarelli et al. (1998); 7- Bacmann 
et al. (2003); 8- Parise et al. (2006); 9- Vastel et al. (2003); 
10- Vastel et al. (2004); 11- Parise et al. (2002); 12- Parise 
etal. (2004); 13- Marcelino et al. (2005). 

there is considerable debate on whether formaldehyde and 
ammonia are formed in the gas phase or on grain surfaces 
(Tielens, 1983; Ceccarelli et al, 2001; Roueff et al, 2005); 
Section 5). From a purely observational point of view, Fig- 
ure 4 points to a similar origin for the three molecules, but 
this is very likely an over-simplification, and maybe just 
due to the common origin of the molecular deuteration, the 
molecular ion H 2 D + (Section 5). 

Indeed, it has been long believed that the deuterium frac- 
tionation is started by the fast ion-neutral reaction H^+HD 
leading to H 2 D+ (Section 5.1). For this reason, H 2 D + has 
long been searched for in regions of massive star formation, 
where, however, it has never been detected (van Dishoeck 
et al., 1992; Pagani et al., 1992). The first detection of 
H 2 D+ was finally obtained by Stark et al. (1999) towards 
the low-mass Class source NGC1333-IRAS4. However, 
the break-through detection is considered to be that ob- 
tained in the Prestellar Core LI 544 (Caselli et al., 2003), 
where the ground ortho-H 2 D + line at 372 GHz is ten times 
brighter than in NGC1333-IRAS4. The observed flux im- 
plies an abundance of H 2 D + similar to that of electrons, 
and, therefore, to that of H^. This was unexpected and un- 
predicted by the chemical models at the time of the dis- 
covery. It is now accepted that the H 2 D + /Hg ratio can 
reach unity in cold gas depleted by heavy-element-bearing 
molecules, as occurs at the centers of Prestellar Cores (Sec- 
tion 2.1). Actually, in the extreme conditions of the cen- 
ters of Prestellar Cores, even the doubly deuterated form 
of Hj^ has been detected (Vastel et al, 2004). The map 
of H 2 D+ emission in L1544 shows that it is very abun- 
dant in a zone with diameter ~ 6000 AU (Vastel et al., 
2006). This zone coincides with the central density plateau, 
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Fig. 4. — The measured N 2 D + /N 2 H + (stars; from Crapsi 
et al. (2005a)), (D 2 CO/H 2 CO) 1/2 (diamonds; Bacmann 
et al. (2003)) and (ND 3 /NH 3 ) 1/3 (triangles; Roueff et al. 

(2005) ) in the sample of Prestellar Cores of Table 2, as func- 
tion of the CO depletion factor. 

where CO is highly depleted (Caselli et al, 1999). This 
further proves the key role played by the freezing of heavy- 
element-bearing molecules in the enhancement of molecu- 
lar deuteration. Finally, one has to note that LI 544 is not a 
peculiarity. An on-going project shows that Prestellar Cores 
are in general strong emitters of the ground ortho-H 2 D + 
line (Tmb > 0.5 K; Caselli et al, in prep.). 

Thanks to the above observations, the very high diag- 
nostic power of the ground-state ortho-H 2 D + line is now 
realized. This is true not only for chemistry, but also for the 
study of the evolutionary and dynamical status of Prestel- 
lar Cores. For example, the ortho-to-para ratio of H 2 D + 
is predicted to critically depend on the density, tempera- 
ture and grain sizes (Flower et al, 2005, 2006). Compar- 
ing these model predictions with observations, Vastel et al. 

(2006) found that the average grain sizes in LI 544 are larger 
than in the ISM (in agreement with the discussion in Section 
2.1). Finally, since very likely most of the heavy-element- 
bearing molecules are frozen out onto the grain mantles in 
the centers of the Prestellar Cores, H 2 D + is the best, if not 
the only way to probe those regions. Of particular inter- 
est is knowing the kinematical structure and whether and 
when the collapse starts at the center, van der Tak et al 
(2005) studied the lines profile of the fundamental ortho- 
H 2 D+ line towards L1544 with this scope. They found that 
the infall likely started at the center of L1544. We predict 
that similar studies with submillimeter interferometers will 
be extremely useful to unveil the first instants of collapse. 

3.2. Class sources 

As for Prestellar Cores, recent observations have re- 
vealed a zoo of molecules with enhanced deuterium frac- 
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Fig. 5.— The measured (D 2 CO/H 2 CO) 1/2 (stars) and 
CH2DOH/CH3OH (squares) in a sample of Class sources 
(Parise et al, 2006) as function of the CO depletion factor. 

tionation. Actually, the discovery of abundant D 2 CO in 
a Class source started the hunt for multiply deuterated 
molecules. Although doubly deuterated formaldehyde had 
been detected fifteen years ago in Orion by Turner (1990), 
the measured abundance (~ 0.2% H2CO) did not seem to 
draw particular attention (probably because Orion has al- 
ways been considered a peculiar region). The discovery that 
D2CO has an abundance about 5% that of H2CO in a low 
mass protostar, IRAS 16293-2422 (Ceccarelli et al, 1998), 
on the contrary, started a flurry of activity in this field, and 
initiated the hunt for multiply deuterated molecules. In a 
few years, the deuterated species quoted in Table 2 were 
detected. 

The best studied Class source is IRAS 16293-2422, 
where multi-frequency observations of singly and doubly 
deuterated formaldehyde constrain the abundance ratios of 
these molecules quite well (Loinard et al., 2000) to be 
HDCO/H2CO = 0.15 and D 2 CO/H 2 CO = 0.05. Besides, 
D 2 CO is very abundant across the entire envelope, which 
extends up to a radius of more than 5000 AU (Ceccarelli 
et al, 2001). Remarkably, the D 2 CO/H 2 CO ratio reaches 
a peak of 0.16 at about 2000 AU from the center. Note 
that the dust at this distance is warm enough to prevent the 
condensation of CO, and the gas is too warm, <~ 35 K, to 
favor the deuterium fractionation in the gas phase (Cecca- 
relli et ai, 2000a; Schoier et ai, 2002). In a few words, 
the two key factors for the enhanced fractionation, cold 
and CO-depleted gas, are absent in the largest fraction of 
the IRAS 16293-2422 envelope, where D 2 CO is abundant. 
Therefore, the observed enhanced deuterium fractionation 
is unlikely to be a present-day product (Ceccarelli et ai, 
2001). In support of this view, the formaldehyde abundance 
profile across the IRAS 16293-2422 envelope increases by 
about a factor of ten where the dust temperature exceeds 50 
K, suggesting that gaseous formaldehyde is due to the par- 
tial, but continuous sublimation of the grain mantles across 



the entire envelope. Layered like onions, the components 
of the mantles sublimate at different distances, correspond- 
ing to different dust temperatures. In the regions where the 
dust temperature is larger than about 25 K, CO-rich mantles 
sublimate, injecting into the gas phase the CO and the trace 
molecules trapped in the ice, like formaldehyde. When the 
dust temperature reaches 50 K the H 2 CO-rich ices subli- 
mate and so on (Ceccarelli et al., 2001). The observed 
deuteration is therefore a heritage of the pre-collapse phase 
(Section 3.1). 

This can also be seen in Figure 5, where the measured 
D/H ratios of formaldehyde and methanol are reported as 
function of CO depletion in a sample of Class sources 
(Parise et al, 2006). Contrary to the Prestellar Cores, there 
is no clear correlation between the molecular deuteration 
and the present CO depletion. The most plausible expla- 
nation is that the deuterated formaldehyde and methanol 
have been built up during the previous phase and that Class 
sources are young enough for it not being substantially 
modified, likely younger than ~ 10 5 yr based on modeling 
(Charnley et al., 1997). Figure 5 also shows that methanol 
is systematically more enriched with deuterium than is 
formaldehyde. This may reflect a different epoch of forma- 
tion of these two molecules on the ices: the formaldehyde 
perhaps formed in an earlier stage than methanol, when CO 
depletion was less severe. Alternatively, present day gas- 
phase reactions may have changed the initial mantle com- 
position already, or, finally, formaldehyde and methanol fol- 
low different routes of formation (gas phase versus grain 
surface?) (Parise etal, 2006). 

Although our comprehension of the mechanisms lead- 
ing to molecular deuteration in the ISM has undoubtedly 
improved in these last few years, there still remain a num- 
ber of unresolved questions. First, the relative role of gas 
phase versus grain surface chemistry is uncertain: it is not 
completely clear yet what molecules and to what extent are 
formed in the gas phase versus on the grain surfaces. Deu- 
terium fractionation promises to be a key aspect in this is- 
sue. Another open question is represented by deuterated 
water. Searches for solid HDO have so far been in vain, 
giving some stringent upper limits to the water deutera- 
tion: < 2% (Dartois et al, 2003; Parise et al, 2003), con- 
siderably less than the deuteration observed in the other 
molecules. One possible explanation is that observations 
of solid H 2 and HDO may be "contaminated" by the con- 
tribution of the molecular cloud in front of the observed 
protostar (see for example the discussion in Boogert et al 
(2002)). A more stringent constraint is set by observations 
of gas-phase HDO in the outer envelope and hot corino, 
where the icy mantles sublimate. Observations towards 
IRAS 16293-2422 set a stringent constraint to the deutera- 
tion of water; in the sublimated ices it is ~ 3%, whereas in 
the cold envelope it is < 0.2% (Parise et al, 2005a). There- 
fore, water deuteration seems to follow a different route 
from formaldehyde, methanol, and the other molecules. 
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4. HOT CORINOS 

While the outer envelopes of Class sources are so cold 
that heavy-element-bearing molecules are at least partially 
frozen onto the grain mantles, the innermost regions are 
characterized by warm dust, so warm that the mantles sub- 
limate totally (Section 2.2). These regions share several as- 
pects of the hot cores, discovered in the 80s around massive 
protostars (Kurtz et al, 2000). Being "smaller" in size and 
mass, they have been called hot corinos (Ceccarelli, 2004b; 
Bottinelli et al, 2004b). Based on the analysis discussed 
in Section 2.2, the radii of the hot corinos are predicted to 
be around 50 AU (IRAS 16293-2422 represents here an ex- 
ception with a radius of 150 AU), comparable to the sizes 
of the Solar System. The key question, and the most rele- 
vant for the chemical structure of Class sources, is which 
molecules are found in the hot corinos, and, in particular, 
whether complex organic molecules are formed, in analogy 
with massive hot cores. More generally, are the hot corinos 
a scaled version of the hot cores, or does a peculiar chem- 
istry take place? And, how reliable are the hot corino sizes 
derived by the analysis of single-dish data? 

Once again, the benchmark for these studies has been 
IRAS 16293-2422. In the hot corino of this source, complex 
organic molecules were detected for the first time (Cazaux 
et al., 2003), and this hot corino was the first to be imaged 
(Kuan et al., 2004; Bottinelli et al, 2004b). In practice, 
all the complex molecules typical of massive hot cores that 
have been searched for were detected, with abundances 
similar to if not larger than those measured in hot cores 
(Cazaux et al, 2003). These species include O and re- 
hearing molecules such as formic acid, HCOOH, acetalde- 
hyde, CH 3 CHO, methyl formate, HCOOCH 3 , dimethyl 
ether, CH3OCH3, acetic acid, CH3COOH, methyl cyanide, 
CH 3 CN, ethyl cyanide, C 2 H 5 CN and propyne, CH 3 CCH. 
Two years after the first detection in IRAS 16293-2422, 
complex organic molecules were also detected towards 
the other Class sources NGC1333 IRAS4A and B, and 



NGC1333 IRAS2 (Bottinelli et al., 2004a; J0rgensen et al, 
2005a; Bottinelli et al, 2006b). The abundances of the 
O-bearing complex molecules detected in the four hot cori- 
nos are shown in Figure 6 as functions of the methanol 
abundance (Bottinelli et al, 2006b). Note that the abun- 
dances have been normalized to the methanol abundance 
so that they do not depend on the uncertain sizes of the 
hot corinos. The number of hot corinos where complex or- 
ganic molecules have been detected is evidently much too 
small to try to draw firm conclusions based on statistical 
considerations. Nevertheless, it is obvious that complex 
organic molecules are rather common in hot corinos. In 
this respect, the prototypical Class source, IRAS 16293- 
2422, may actually be representative for Class sources. 
Furthermore, the abundances of HCOOH, HCOOCH3 and 
CH3OCH3 in hot corinos are comparable to the methanol 
abundance and to formaldehyde (Bottinelli et al, 2006b), 
and do not depend appreciably on the source luminosity, 
nor on the methanol and formaldehyde abundance. If the 
O-bearing complex molecules are formed in the gas phase 
from methanol and formaldehyde, as some theories pre- 
dict (Section 5.3), then they "burn" the majority of their 
"parent" molecules, and this burned fraction is independent 
of the amount of H 2 CO and/or CH3OH sublimated from 
the grain mantles. Alternatively, if the O-bearing complex 
molecules are synthesized on the grain surfaces, they are an 
important ice component. So far, however, only HCOOH 
has been claimed to be found in its solid form (Schutte 
et al, 1999; Gibb et al, 2000). The other evident message 
of Figure 6 is that hot corinos are not a simply scaled ver- 
sion of hot cores. The abundance of the O-bearing complex 
molecules in hot corinos is larger than that in hot cores, by 
more than a factor of ten: HCOOH is even 100 times more 
abundant. 

So far, we have focused on results based on single-dish 
observations. These observations encompass regions of at 
least 1000 AU, whereas the predicted sizes of the hot cori- 
nos do not exceed 150 AU. Recently, interferometric ob- 
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servations have resolved the emission of complex organic 
molecules in the hot corinos towards IRAS 16293-2422 
(Kuan et al, 2004; Bottinelli et al, 2004b; Chandler et al., 
2005), NGC1333-IRAS2 (J0rgensen et al, 2005a), and 
NGC1333-IRAS4B (Bottinelli et al. , in prep.). From these 
data, a few general conclusions can be drawn. First, for 
IRAS 16293-2422 and NGC1333-IRAS4B, the line emis- 
sion from the complex molecules is concentrated in two 
spots, centered on the two objects forming in each case a 
proto-binary system (Wootten, 1989; Looney et al, 2003; 
Sandell and Knee, 2001; Di Francesco et al., 2001). The 
full line flux of the observed complex molecules, detected 
with the single-dish, is recovered in these two spots, which 
means that all the single dish emission originates in those 
regions. In addition, there is no sign of outflows in the line 
emission of complex molecules; the emission is compact 
and the line profile is relatively narrow. The brightest - and 
so far only resolved spot - in IRAS 16293-2422 has a radius 
of ~ 150 AU, close to the size expected for warm (T > 100 
K) gas (Ceccarelli et al., 2000a). As emphasized above, 
each of these two objects is indeed a protobinary system. 
In each case, the two objects in the protobinary have differ- 
ent characteristics; eg., one spot is bright in the continuum 
but barely detected in the complex molecule line emission, 
while this is reversed for the other spot. This difference 
does not necessarily imply that the abundances in the indi- 
vidual spots of these protobinaries are truly different. The 
spot with the weaker line emission is also unresolved and 
the line likely optically thick (Bottinelli et al, 2004b). Ev- 
idently, high-resolution observations of complex organic 
molecules are a valuable tool to study the chemistry of the 
hot corinos, and how it depends on physical conditions, 
rather than the mantle composition and age, which presum- 
ably are the same for the two objects of the binary system. 

5. CHEMICAL MODELS 

5.1. Molecular deuteration in the gas phase 

Although the basic gas-phase mechanism of molecular 
deuteration in the cold interstellar medium was elucidated 
almost thirty years ago (Guelin et al, 1977; Watson et al, 
1978) and incorporated into detailed models over a decade 
ago (Millar et al, 1989), it is only recently that an un- 
derstanding of the effect in star-forming regions has been 
achieved by astrochemists. Since the reservoir of deuterium 
lies in the molecule HD, deuteration proceeds by removing 
deuterium atoms from this molecule onto the many other 
trace constituents of the gas. The principal mode of trans- 
fer at low temperatures has been thought to be the reaction 
system 

H++HD^H 2 D++H 2 , (2) 

in which the left-to-right reaction is exothermic by approx- 
imately 230 K. If no other reactions need be considered to 
determine the abundance ratio of H2D + to Hj", the slowness 
of the backwards reaction leads to the stunning prediction 
that the abundance of the deuterated ion, now referred to as 



an isotopologue, can exceed that of Hg . For cold cores at 
the standard density of 10 4 cm~ 3 , the prediction is in error 
because H 2 D + reacts quickly with both electrons and heavy 
species such as CO. The reactions with heavy species serve 
to spread the enhanced deuteration around so that typical ra- 
tios between singly-deuterated isotopologues (e.g., DCO + ) 
and normal species (e.g. HCO + ) are predicted to lie in the 
range 0.01-0.10, in agreement with observation. Two other 
exchange reactions are of lesser importance: these involve 
the ions CH^ and C2H;}" . Their deuterium exchange reac- 
tions with HD leading to the ions CH2D+ and C2HD+ are 
more exothermic than the H3 case, so that the backwards 
reactions are slow up to higher temperatures, and subse- 
quent reactions involving the deuterated ions can maintain 
deuterium fractionation when H 2 D+ is no longer enhanced. 
Other reactions are also non-negligible, including some in- 
volving deuterium atoms (Brown and Millar, 1989). 

As the density increases, the picture changes dramati- 
cally. Near the center of a cold pre-stellar core, for exam- 
ple, the gas density increases to 10 5 ~ 6 cm~ 3 and the time 
scale for accretion of species onto dust particles becomes 
so short that most heavy species are strongly depleted onto 
dust particles. In addition, the fractional ionization in the 
gas is thought to become quite low if one takes account of 
gas-grain interactions to deplete the chemically rather inert 
atomic ion H + (Liszt, 2003). Since seemingly all destruc- 
tion mechanisms of H 2 D+ are slow, a large abundance of 
H 2 D+ does indeed build up, in agreement with observations 
(Section 3.1). But the story does not end here: in analogy 
to the reaction above, further reactions with HD convert a 
significant portion of the H 2 D+ to the more deuterated iso- 
topologues HD^ and even D^ (Roberts et al, 2003; Walm- 
sley et al, 2004; Ceccarelli and Dominik, 2005). Indeed, 
simple chemical models of the densest portions show the 
triply deuterated ion to be the most abundant of the four 
ions in the series, and possibly the most abundant ion in the 
gas (Roberts et al, 2003). Figure 7 shows the abundances 
of the deuterated forms of H^", H + and e~, as functions of 
the CO depletion factor respectively. Since the age and any 
other parameter is constant, increasing the depletion factor 
corresponds to increasing the density. For "intermediate" 
(< 30) CO depletion factors, H 2 D+ is the most abundant 
deuterated form of H^. However, for larger depletions D^ 
takes over, and becomes the dominant positive charge car- 
rier. A similar effect leads to the ions CD^ and C2D2^ . 

Although a tiny region in the center may be completely 
devoid of heavy species (Walmsley et al, 2004), such 
molecules do exist, albeit with depleted abundances, near 
the center (Bergin et al, 2002; Pagani et al, 2005; Lis 
et al, 2006). The remaining heavy species will react with 
all three of the deuterating ions, especially D^, leading to 
highly-deuterated isotopologues such as D 2 CO and ND 3 . 
Critical here is the high efficiency of D J since, unlike the 
partially deuterated analogs, it can only deuterate and not 
protonate heavy species. Contrast the reaction 

D+ + CO — > DCO+ + D 2 (3) 
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with 

H 2 D+ + CO — >DCO + + H 2 ;HCO++HD. (4) 

In the latter reaction, deuteration occurs on only one of three 
collisions. The increase in efficiency with Dg as the deuter- 
ating agent is especially important in the formation of mul- 
tiply deuterated neutral isotopologues, which occurs via a 
cycle of deuteration followed by dissociative recombina- 
tion reactions. For example, starting with normal ammo- 
nia, deuteration leads to NH3D+, followed by dissociative 
recombination to form NH 2 D, followed by deuteration to 
form NH 2 D 2 f , followed by dissociative recombination to 
form NHD 2 , etc (Roueffet al, 2005). 

Gas-phase chemical models of Prestellar Cores includ- 
ing deuteration range from simple homogeneous treatments 
relevant to the center of the core {Roberts et al, 2003) to 
shell models that can be static (Roberts et al, 2004) or in a 
state of collapse (Aikawa et al, 2003, 2005). Accretion onto 
grains must be included in all approaches. The collapse can 
be treated via analytical approaches (Aikawa et al., 2003) 
or hydrodynamic simulations (Aikawa et al, 2005). The 
chemistry occurring on grain surfaces can also be included 
(Aikawa et al, 2003, 2005). A relatively simple static ap- 
proach by Roberts et al. (2004) illustrates the power of the 
chemical models. In terms of column densities, Roberts 
et al. (2004) obtain a value for D 2 CO higher by only a fac- 
tor of three than that observed by Bacmann et al. (2003), 
and a CO depletion of a factor of 20 compared with the 
standard value, in good agreement with the observed deple- 
tion of 14 (Bacmann et al, 2002) (see also Section 3 and 
2). The overall computed D 2 CO/H 2 CO ratio is 0.3, also in 
good agreement with observation. In general, models with 
collapse are not in as good agreement with observation be- 
cause it is not facile to vary the details of the collapse to 
obtain optimal agreement. 



5.2. Mantle formation and surface chemistry 

As the gas-phase chemistry occurs and mantles accrete 
in cold sources, a surface chemistry is also active. A proper 
mathematical treatment of diffusive surface reactions on in- 
terstellar grains is not simple. The first such treatments used 
rate equations analogous to those used for gas-phase pro- 
cesses (Pickles and Williams, 1977), and the method is still 
the only practical approach to the construction of large gas- 
grain models (Ruffle and Herbst, 2000). The method is only 
a rough approximation to reality for a number of reasons, 
and a variety of more accurate but more complex methods 
have been advocated. If we imagine grains to be rather 
smooth surfaces on which binding parameters for adsor- 
bates do not differ strongly from one site to another, there 
is the problem that only small (fractional) numbers of reac- 
tive species such as atomic hydrogen are present on a given 
grain on average. With such small numbers, it is more rea- 
sonable to consider discrete rather than average quantities. 
Both the discreteness of the problem and the large fluctua- 
tions possible from grain to grain then argue for a stochastic 
approach. Following the pioneering work of Tielens and 
Hagen (1982), a truly stochastic method was utilized by 
Charnley (2001) based on Monte Carlo methods, while a 
master equation approach was advocated by Biham et al. 
(2001) and by Green et al (2001). The latter approach has 
recently been incorporated into moderately sized gas-grain 
networks by Biham and Lipshtat (2002), Lipshtat and Bi- 
ham (2004) and Stantcheva and Herbst (2004). 

To add to the complexity of the problem, interstellar sur- 
faces are likely to be irregular, porous, and even amorphous, 
so that binding sites for adsorbates are likely to be diverse, 
including some that contain much deeper wells than others. 
A new extension of the Monte Carlo method for irregular 
grains with silicate, carbonaceous and icy surfaces has been 
introduced by Chang et al (2005) and Cuppen and Herbst 
(2005) for the specific case of the formation of hydrogen 
molecules. Here, the fraction of binding sites with deep 
potential wells totally changes the character of the prob- 
lem since numerous H atoms are present over wide tem- 
perature ranges on each grain. Whether or not an approach 
of this complexity can be extended to large models is un- 
clear. A simpler approach to different binding sites is the 
physisorption-chemisorption model of Cazaux and Tielens 
(2004), where only two types of binding sites are consid- 
ered for H atoms, shallow and deep. In this approximation, 
rate equations can be used to a reasonable degree of accu- 
racy. 

At later stages, the surface chemistry includes deuter- 
ation via atomic deuterium, which is produced in the gas 
by dissociative recombination of deuterated ions; viz., 
+ c~ — ► D + D + D, and can lead in very dense re- 
gions to an abundance of atomic deuterium nearly equal 
to the abundance of atomic hydrogen (Roberts et al, 
2003, 2004). The H and D atoms land on grains and, 
because of their rapid diffusion and high surface re- 
activity, both hydrogenate and deuterate heavy atoms 
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and reactive molecules to form more saturated forms 
via sequential reactions. For example, the formation of 
water ice can occur by the addition of two hydrogen 
atoms to an oxygen atom that lands on a grain surface: 
O + H — ► OH; OH + H — > H 2 0. or can start from OH, 
which can be formed on grain surfaces by successive re- 
actions involving O2 and O3 (Helens and Hagen, 1982). 
Note that even if the gaseous O2 is low, an appreciable 
abundance of O2 can built up on the surface through reac- 
tions of O plus O, or 2 can be formed in the reaction of 
O3 with H. There is even laboratory evidence that hydro- 
genation of CO into formaldehyde (H 2 CO) and methanol 
(CH3OH) occurs (Watanabe et al, 2003; Hidaka et al, 
2004): CO -> HCO H 2 CO -> H 3 CO CH 3 OH , al- 
though two of the hydrogenation reactions possess small 
activation energy barriers. Following a similar path, not 
only H 2 CO and CH3OH are formed but also their deuter- 
ated isotopologues (Charnley et al, 1997). It is interesting 
to note that there is currently no viable formation mecha- 
nism for methanol in the gas phase, so that all methanol 
detected in the gas has its genesis on grain surfaces. This 
statement is true even for the small amount of methanol 
detected in cold sources. Here surface formation is fol- 
lowed by inefficient desorption. The results of a number 
of models - which differ in the sophistication of the gas- 
grain interaction, their assumptions on initial conditions 
including gas-phase abundances, and the chemical routes 
included - have been reported in the literature (Helens and 
Hagen, 1982; Hasegawa et al., 1992; Aikawa et al., 2003; 
Stantcheva and Herbst, 2003; Aikawa et al, 2005). Deuter- 
ation on cold surfaces was first explored by Helens (1983) 
and later studied by Brown and Millar (1989), Charnley 
et al. (1997), Caselli et al. (2002a), and Stantcheva and 
Herbst (2003) among others. The basic picture adopted by 
astrochemists is that D atoms can deuterate atoms and radi- 
cals in a similar manner to the reactions of H atoms. If we, 
for example, consider the competing deuteration and hy- 
drogenation processes starting with surface CO, all of the 
different isotopologues and isotopomers of methanol can 
be produced: CH 3 OD, CH 2 DOH, CH 2 DOD, CHD 2 OH, 
CHD 2 OD, CD3OH, and CD 3 OD (Charnley et al, 1997). 
Their relative abundances on grain surfaces depend on the 
flux ratio between H and D atoms landing on the grains. 
That the deuteration is more complex than envisioned by 
astrochemists, however, is shown by some recent experi- 
mental evidence provided by Nagaoka et al. (2005) where 
it seems that even normal methanol can be deuterated by 
reactions with atomic D on grain surfaces. 

All model calculations confirm that ices are produced by 
surface reactions, and that these ices are, at least on aver- 
age, dominated by water, as observed. The large amount of 
CO ice detected by observers is probably mainly the result 
of accretion from the gas, where this molecule is produced 
copiously. The production of surface methanol occurs once 
there is sufficient CO to be hydrogenated. The production 
of the other very abundant surface species - carbon dioxide 
(C0 2 ) - is not well understood. 



5.3. Mantle evaporation and complex organic molecule 
formation 

A variety of processes can return ice mantle species to 
the gas phase, even during the cold prestellar core phase. 
The most important one here is cosmic ray driven desorp- 
tion where an Fe-cosmic-ray (E ~ 10 — 100 MeV/nucleon) 
temporarily heats the grain to a high enough temperature 
that some mantle species evaporate (Leger, 1983; Leger 
et al, 1985). By itself, this process will only return weakly 
bound species (eg., CO, N 2 ) to the gas phase. However, 
this process becomes very effective, even for more tightly 
bound species such as H 2 0, if the ices contain chemical 
energy in the form of stored radicals produced by UV pho- 
tolysis (D'Hendecourt et al, 1982, 1986; Helens, 2005). 
It should be emphasized that the thermal shock driven into 
the grain by the penetrating cosmic ray can also "lift" off a 
small quantity of larger species (Johnson et al, 1991). 

As the protostellar stage commences and temperatures 
begin to increase, the "dirty" grain mantles will be lost par- 
tially or totally via evaporation in a rather complex process 
stemming from the heterogeneous nature of the mantles 
(Viti et al, 2004). Nearer to the bipolar outflows associ- 
ated with protostellar regions, shock waves can also lead 
to the removal of the mantles. In particular, heating of ice 
mantles will lead to evaporation and an injection into the 
warm dense gas of large quantities of grain surface pro- 
duced species such as water, methanol, and formaldehyde. 
Reactions among these simple molecules can then lead to 
more complex species (Millar et al, 1991; Charnley et al, 
1992; Caselli et al, 1993; Rodgers and Charnley, 2003; 
Nomura and Millar, 2004). The process begins with the 
production of H^ via cosmic ray bombardment of H 2 , fol- 
lowed by the reaction of H^ and H 2 . This triatomic species, 
possibly through some intermediaries, transfers its proton 
to methanol. It has long been believed that the resulting 
protonated methanol, CHsOH^ , will rapidly undergo alkyl 
cation transfer reactions with other species. With CH3OH 
and H 2 CO, this might lead to the formation of dimethyl 
ether and methyl formate, which are abundant complex 
molecules in hot corinos (Section 4). However, the above 
reactions advocated to produce methyl formate do not oc- 
cur, based on laboratory experiments (Horn et al, 2004). 
A more successful synthesis seems to be the one leading to 
dimethyl ether: 

CH3OH + CH3OH+ — ► CH3OHCH+ + H 2 0, (5) 

CH3OHCH+ + e" — ► CH3OCH3 + H, (6) 

but it is to be noted that experimental measurements of 
the neutral products of dissociative recombination reactions 
show that two-body products such as shown here are often 
minor channels. Eventually, these complex species are de- 
stroyed on a timescale of 10 4-5 yr, depending on the tem- 
perature and density of the gas, to reform CO. Indeed, if 
the chemistry were understood well enough, abundances of 
these species have the potential to be a chemical "clock", 
timing the formation of the central object, which has heated 



11 



up and evaporated the ices. Other determinants of age are 
abundance ratios of sulfur-bearing species, although the use 
of these indicators requires sophisticated modeling of both 
the chemical and the physical structure of the studied ob- 
ject (Wakelam et ai, 2004, 2005). The need for such a so- 
phisticated approach arises because the abundance ratios of 
sulfur-bearing species, used as clocks, turn out to also be 
dependent on initial abundances following evaporation as 
well as on physical conditions. 

Deuterium fractionation in the ices arises mainly from 
reactions involving deuterium atoms accreted from the 
highly fractionated gas. Such fractionation occurs in the 
abundant surface species (H 2 0, H 2 CO, CH 3 OH, NH 3 , 
CH4) acting as signposts of a cold prestellar-core history in 
which gaseous D atoms are abundant. When these species 
are returned to the gas, the fractionation becomes more 
apparent. Although it is sometimes difficult to disentan- 
gle purely gas-phase fractionation from surface fraction- 
ation followed by evaporation, fractionation in gaseous 
molecules formed mainly on surfaces (e.g., methanol) is 
due at least initially to surface reactions. This fractiona- 
tion will also be passed on to daughter products such as 
dimethyl ether and methyl formate via gas-phase processes. 
Thus, the deuteratium fractionation can be used as a tracer 
of the chemical routes involved in the formation of complex 
molecules in hot cores. Indeed, the fractionation pattern of 
dimethyl ether and methyl formate should reflect directly 
that of methanol and formaldehyde. Eventually, reactions in 
the warm dense gas of the hot core will reset this deuterium 
fractionation to that appropriate for the gas temperature 
(Charnley et ai, 1997). This chemical kinetic evolution 
may run differently for different species or even for differ- 
ent chemical groups in one species. Specifically, for the 
fractionation pattern of methanol, it has been suggested 
that protonation of methanol isotopomers followed by dis- 
sociative electron recombination reforming methanol will 
preferentially remove the deuteration from the OH group 
but not affect the CH3 group {Osamura et ai, 2004). The 
analysis may be incorrect, however, because the most re- 
cent storage ring experiments (Geppert et ai, 2004) show 
that protonated and deuterated methanol ions undergo dis- 
sociative recombination reactions that lead to virtually no 
neutral methanol or its isotopologues at all. If all of the gas- 
phase isotopologues of methanol are depleted at the same 
rate, then the abundance ratios are indeed determined by 
the surface chemistry during the previous cooler eon. 

6. CONCLUSIONS: CLOSED AND OPEN QUES- 
TIONS 

In this contribution, we have summarized new insights 
into the physical and chemical composition and evolution 
of Prestellar Cores and Class low-mass Protostars. In par- 
ticular, 

1. On the physical structure, Prestellar Cores can be 
described by the density distribution of a Bonnor-Ebert 
sphere, whereas Class sources are consistent with the Shu 



"inside-out" recipe. 

2. At the center of Prestellar Core condensations, den- 
sities are so high and temperatures so low that possibly all 
heavy-element-bearing molecules freeze out onto the dust 
grains, forming mantles of ices. We discussed the possibil- 
ity that small grains coagulate into larger grains during this 
phase. 

3. A similar molecular depletion is observed in large 
regions of the outer envelopes of Class sources. These re- 
gions are indistinguishable from Prestellar Cores on this ba- 
sis, arguing that the latter are precursors of Class sources. 

4. The molecular depletion is accompanied by a dra- 
matic enhancement of the molecular deuteration, observed 
in both Prestellar Cores and Class sources. In the in- 
nermost regions of Prestellar Cores, the only observable 
molecules may be H 2 D+ and HD^, which are likely as 
abundant as Hg" or even more so. 

5. In Class sources, many multiply deuterated 
molecules have been observed. Even triply deuterated 
molecules have been detected, with abundances enhanced 
by up to 13 orders of magnitude with respect to the D/H 
elemental abundance. 

6. Triggered by the observations of multiply deuterated 
molecules in Prestellar Cores and Class sources, a new 
class of models for deuteration has been developed in the 
last few years. These models predict that the low temper- 
atures coupled with the disappearance of heavy-element- 
bearing molecules from the gas phase cause the forma- 
tion of abundant H 2 D + , HD^" and D^, which are ex- 
tremely efficient in passing their deuterium atoms on to 
other molecules. In extreme conditions, is predicted to 
be the most abundant molecular ion in the gas, surpassing 
by orders of magnitude the number density of H3 . 

7. During the cold phases of low-mass star formation, 
grain mantles are formed, consisting mainly of water ice, 
but with many other molecules in quantities that can be im- 
portant. Extreme molecular deuteration is a clear hallmark 
of this phase and is recorded in the mantles. Hydrogenated 
molecules, like formaldehyde and methanol, are also be- 
lieved to be formed on the grain surfaces during the cold 
pre-collapse phase. 

8. Class envelopes consist of two chemically dis- 
tinct regions: the outer envelopes and the inner regions, 
called hot corinos. The border between the two is where 
the dust temperature reaches the sublimation temperature 
of the grain ices. This has a major impact on the gas chem- 
istry. In the hot corinos, the chemistry is dominated by the 
evaporation of the mantles, built up during the pre-collapse 
phase. Once the mantle components are in the gas phase, 
they undergo successive reactions leading to the formation 
of many complex organic molecules, similar to those ob- 
served in the hot cores. 

Although important progress has been achieved in these 
last few years, it has also been accompanied by the rise 
of numerous new questions. For example, while molecu- 
lar depletion is now an accepted "fact", it is still not to- 
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tally clear whether in the innermost regions of the Prestellar 
Cores the density increases so much that not even a trace 
of heavy-element species remains in the gas phase. Also, 
it still remains unclear when and exactly how the collapse 
starts. It is clear now that observations of heavy-element- 
bearing molecules may not be the best way to proceed. In 
contrast, as we now realize, the ground-state transition of 
ortho-H 2 D + (and, to a lesser extent, para-HD^ ) provides a 
new and exciting probe with which to search for the long- 
sought collapse in the innermost regions of Prestellar Cores 
(van der Tak et al, 2005). However, in order to fully ex- 
ploit this new diagnostic, we have to first fully understand 
the ortho-para ratio of H 2 D + (Flower et al, 2006), which 
is not measured at all so far. It is also clear that mod- 
els for the Prestellar Core phase have improved dramati- 
cally with the inclusion of the extreme depletions of the gas 
phase and the formation of the multiply deuterated forms of 
Hjj~. However, these models still struggle to reproduce the 
extreme observed deuteration, especially that of formalde- 
hyde and methanol. While this is likely due to grain surface 
chemistry, it remains the case that, unless there is an en- 
hancement in the reactivity of D over H, the gas landing 
on the grain surfaces at the moment of the mantle forma- 
tion must have an atomic D/H ratio higher than accounted 
for by the current gas phase chemistry models. We also 
have little knowledge as to whether the prediction that D^ 
is the most abundant molecular ion in extreme conditions 
is correct. The coming years will likely see many of these 
questions answered, with the advent of the satellite Her- 
schel and, foremost, the interferometer ALMA. 

If many questions remain about the cold pre-collapse 
phase, even more questions are unanswered for the proto- 
stellar phase, and many have likely not even been asked. 
For example, the very nature of the hot corinos is still very 
much debated. Several authors question how much of the 
observed warm gas, with complex organic molecules, is 
due to the interaction of an outflow with the inner envelope 
(Chandler et ah, 2005). Another possibility invoked is that 
this warm gas enriched with complex organic molecules re- 
sides in a disk-like atmosphere (J0rgensen et al., 2005a). 
The reality is indeed that very few observations have re- 
solved hot corinos; actually, in just one source has this so far 
occurred (Kuan et al, 2004; Bottinelli et al, 2004b). Un- 
til more sensitive observations are available with SMA and 
ALMA, the question will remain open for debate. Another 
open question regards the chemistry leading to the observed 
molecular complexity in hot corinos. No clear observations 
are able so far to distinguish the roles of gas phase and grain 
surface chemistry in the formation of formic acid or methyl 
formate, just to mention two molecules frequently detected. 
The answer to this question will likely need not only more 
and more sensitive observations towards the hot corinos, but 
also laboratory experiments coupled with modeling in or- 
der to elucidate the role of the grains in the story. Last, but 
not least, we suspect that we are just viewing the tip of the 
iceberg: what is still hidden from our eyes? What is the ul- 
timate molecular complexity reached during the hot corino 



phase? Do biotic molecules form? Do they play a role in 
the trigger and/or diffusion of life in the Universe (see the 
chapter by Gaidos & Selsis)? 

Finally, many more questions exist regarding the fate of 
the molecules formed through the pre-stellar and protostel- 
lar phases. It is likely that some condense again onto the 
grain mantles during the phase of the proto-planetary disk, 
at least in the outer zones of the disk that are shielded by 
FUV photons (see the chapters by Bergin et al, Dutrey et 
al. and Dullemond et al). The same grains are likely to 
coagulate into larger aggregates, which are the bricks from 
which the planets eventually form (see the chapters by Do- 
minik et al. and Wadhwa et al). How much of the mantle 
formed in the Prestellar Core phase is preserved? Do the 
molecules, and particularly the complex organic molecules 
formed during the hot corino phase, freeze-out onto the new 
mantles without other alterations? Or what kind of alter- 
ations do they undergo? How pristine is the material form- 
ing the comets, asteroids and planets? 

To reconstruct the full story, we need many more pieces 
of the puzzle. One strong motivation for pursuing the story 
of protoplanetary disks is that their aftermath may have 
been part of the history of our own Solar System. The water 
of the oceans is an emblematic example. The ratio between 
HDO and H 2 (known as SMOW) is - 10~ 4 (De Witt 
et ah, 1980), fully ten times larger than the D/H elemental 
abundance of the Solar Nebula (Geiss, 1993). The origin of 
the oceans has been long debated and new theories keep 
coming. The deuterium enhancement is in line with the 
theory that much of the ocean water comes from the bom- 
bardment of the early Earth by comets (Owen and Bar-Nun, 
1995), for comets show a similarly enhanced HDO/H2O 
ratio (Meier et al, 1998; Bockelee-Morvan et al, 1998). 
Other possibilities include water brought by smaller bodies 
of the Solar System (meteorites and asteroids) (Raymond 
et al, 2004; Gomes et al, 2005). These bodies, as well 
as comets, have the imprint of the pre-collapse/protostellar 
phase in their ices (see the chapters by Jewitt et al, Wooden 
etal. and Alexander et al). Regarding deuterium fractiona- 
tion, the question remains as to what extent the enhanced 
HDO/H2O ratio is a legacy from the protostellar/proto- 
planetary phase and to what extent it has been modified 
by subsequent chemistry. And in addition to water, what 
other molecules were inherited by the early Earth and how 
did they arrive? As with other molecular indicators, how 
much of the terrestrial HDO/H 2 ratio is a legacy depends 
on how it evolved during the various phases of the earth's 
formation. So far, we have constraints on the evolution of 
deuterium fractionation from the pre-collapse, protostellar, 
and proto-planetary phases for only two molecules, HCO + 
and H 2 0, and in the outer zones of the disk only, which 
may not necessarily be connected with cometary forma- 
tion. The case of DCO + /HCO + , representing deuteration 
in the gas phase, seems to keep a rather constant value, 
around 10% (van Dishoeck et al, 2003). The HDO/H 2 
ratio, very likely a product of the grain mantles, seems to 
remain quite similar too, at a level of ~1%, with just a rel- 
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atively small decrease in the proto-planetary phase (Cec- 
carelli et al, 2005; Dominik et al, 2005). Evidently, this 
information alone is far too little to resolve the problem of 
the oceanic origin, but shows how important the study of 
chemistry in the prestellar and protostellar phases of star 
formation is and how much it will eventually aid our un- 
derstanding of later phases up through the present eon. The 
coming decade will likely see many of these questions an- 
swered with the advent of the satellite Herschel and the sub- 
millimeter interferometer ALMA. 
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